Entanglement in S'0(3)-invariant bipartite quantum systems 
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The structure of the state spaces of bipartite N <g> N quantum systems which are invariant under 
product representations of the group SO (3) of three-dimensional proper rotations is analyzed. The 
subsystems represent partides of arbitrary spin j which transform according to an irreducible rep- 
resentation of the rotation group. A positive map # is introduced which describes the time reversal 
symmetry of the local states and which is unitarily equivalent to the transposition of matrices. It is 
shown that the partial time reversal transformation $2 = / <S> $ acting on the composite system can 
be expressed in terms of the invariant 6-j symbols introduced by Wigner into the quantum theory 
of angular momentum. This fact enables a complete geometrical construction of the manifold of 
states with positive partial transposition and of the sets of separable and entangled states of 4 <g) 4 
systems. The separable states are shown to form a three-dimensional prism and a three-dimensional 
manifold of bound entangled states is identified. A positive maps is obtained which yields, together 
with the time reversal, a necessary and sufficient condition for the separability of states of 4 ® 4 
systems. The relations to the reduction criterion and to the recently proposed cross norm criterion 
for separability are discussed. 

PACS numbers: 03.67.-a,03.65.Ud,03.65.Yz 
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I. INTRODUCTION 

It is one of the bàsic postulates of quantum mechanics 
that the Hilbert space of states of a composite system 
is given by the tensor product of the Hilbert spaces per- 
taining to its subsystems. If a system is composed of 
two JV-state systems with Hilbert space H = C N , the 
mixed states of the total system are represented by den- 
sity matrices which act on the tensor product C N <g> C N . 
A state of such an N ® N system is said to be separable 
or classically correlated if it can be generated by mixing 
with certain probabilities an ensemble of states which 
describe statistically independent subsystems States 
which cannot be represented in this way are called insep- 
arable or entangled. The characterization, classification 
and quantification of mixed-state entanglement and the 
development of explicit necessary and sufficient separa- 
bility criteria turn out to be an extremely hard problem 
■ The solution of this problem could have far-reaching 
consequences for fundamental qüestions of quantum me- 
chanics and computational complexity theory 0] and for 
applications in the theory of quantum information 0. 

A great simplification of the entanglement problem is 
obtained though the introduction of symmetries Q^Q!- 
By the requirement of invariance under certain groups 
of symmetry transformations one restricts the set of all 
states to a low-dimensional manifold of invariant states 
and one may hope to get a tractable problem which 
is solvable with the help of group theoretical and alge- 
braic methods. A prominent example is gjyen by the 
one-parameter family of the Werner states [ï| which re- 
sults from the requirement of invariance under all prod- 
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uct transformations of the form U ® U, where U varies 
over the full group of unitary N x N matrices. A further 
related example is the one-parameter family of isotropic 
states @, Ü| which are invariant under all product uni- 
tàries U ® U* , where U* denotes the complex conju- 
gation of U. Imposing the invariance under all trans- 
formations of the form O (8 O, where O belongs to the 
group of orthogonal NxN matrices, one obtains the two- 
dimensional manifold of orthogonally invariant states 0|. 
It is clear that the larger the symmetry group the smaller 
is the remaining space of invariant states and the easier 
should be the analysis of its structure. In fact, the prob- 
lem of the explicit determination of the separable states 
under symmetry requirements can be solved completely 
for the examples given above. 

A physically natural symmetry group is the group 
50(3) of proper rotations in three dimensions. The un- 
derlying assumption is that the states of the subsystems 
transform according to an JV = (2j + l)-dimensional ir- 
reducible representation of the rotation group which cor- 
responds to a fixed angular momentum j. The subsys- 
tems thus behave under rotations as partides with a cer- 
tain spin j, The rotation group then operates through 
a reducible product representation on the states of the 
composite system. Any 50(3)-invariant state can be 
decomposed into a sum of projections onto the irre- 
ducible subspaces belonging to the various eigenvalues 
J = 0, 1, . . . , 2j of the total angular momentum oper- 
ator. This shows that the rotationally invariant states 
form an (JV — l)-dimensional manifold. The requirement 
of 50(3) invariance reduces the full (A^ 4 — l)-dimensional 
space of all mixed states of an N<g>N system to an (N— 1)- 
dimensional space of invariant states. 

The invariance under 50(3) transformations repre- 
sents in general a much smaller symmetry than those of 
the examples given above. For example, the manifolds of 
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the Werner states and of the isotropic states can be em- 
bedded into the set of rotationally invariant states. These 
examples are thus special cases of the SO(3) symmetry. 

The problem of mixed-state entanglement in £0(3)- 
invariant bipartite systems will be analyzed in this paper. 
We find that the state spaces exhibit an interesting con- 
vex structure and several important phenomena as the 
emergence of non-decomposable positive maps 0] and 
bound entanglement [ïjj- The physical significance of 
the SO(3) symmetry derives from the fact that any rota- 
tionally invariant state can be produced from the maxi- 
mally entangled angular momentum singlet state J = 
through the application of an isotropic dynamical map 
which operates locally on the subsystems. The set of 
S'0(3)-invariant states is thus identical to the set of states 
which results from interactions of the singlet state with 
noisy isotropic environments. 

A powerful tool in studies of entanglement is the op- 
eration of taking the partial transposition of states. The 
requirement of positive partial transposition (PPT) rep- 
resents a strong necessary condition for the separability 
of states, known as the Peres criterion 0|. A conceptu- 
ally simple but crucial point of the present investigation 
consists in the replacement of the transposition by an- 
other unitarily equivalent operation which is identical to 
the time reversal transformation of partides with spin j. 
It is known that the anti-unitary operation of the time re- 
versal commutes with the representations of the rotation 
group. This fact allows one to characterize the partial 
transposition by means of invariant quantities which are 
directly related to Wigner's 6-j symbols [ïïij . 

The 6-j symbols arise in the transformation between 
different coupling schemes for the addition of angular mo- 
menta [Ï3 |. They can be expressed as invariant sums 
over products of vector-coupling (Clebsch-Gordan) coef- 
ficients. Thus we find a close connection between the par- 
tial transposition, the time reversal symmetry and cer- 
tain group-theoretical invariants built out of the vector- 
coupling coefficients. It will be demonstrated here that 
this connection to group-theoretical concepts leads to im- 
portant implications on the entanglement structure of the 
state space. 

The content of the paper can be summarized as fol- 
lows. In Sec.|ïï]we briefly recali some facts from the rep- 
resentation theory of the rotation group and introduce 
an appropriate parametrization of the set of rotationally 
invariant states. The partial transposition and the corre- 
sponding transformation of the partial time reversal are 
discussed in Sec. IIIII It is shown that this transforma- 
tion preserves the rotational invariance of states and its 
relation to the Wigner 6-j symbols is derived. 

These results are used in Secs. IIVI and Ivl to develop 
a geomètric representation of the sets of the PPT states 
and of the separable states in the cases N — 2, 3 and 4. 
Most importantly, in the case N — 4 we find that the set 
of separable states is isomorphic to a three-dimensional 
prism, i. e., to a polyhedron which is bounded by three 
squares and two triangles. We further identify a three- 



dimensional manifold of bound entangled states with pos- 
itive partial transposition. 

Finally, Sec. IVII contains a discussion of the results 
and a number of conclusions which can be drawn from 
the present investigation. In particular, we construct a 
positive map which yields, together with the time rever- 
sal, a necessary and sufficient condition for separability 
in the case of 4 ® 4 systems. Moreover, we discuss the 
relations to two further criteria of separability, namely 
the reduction criterion and the cross norm criterion. 



II. THE SET OF SÜ(3)-INVARIANT STATES 
A. Representations of the rotation group 

We consider a bipartite quantum system whose lo- 
cal parts are iV-state systems with corresponding state 
space H = C N . The Hilbert space of the composite sys- 
tem is given by the tensor product space Ji ®7i. The 
local state spaces are regarded as angular momentum 
manifolds corresponding to a certain eigenvalue of the 
square of the angular momentum operator j. Thus, the 
state space H is spanned by a fixed orthonormal ba- 
sis of N = (2j + 1) angular momentum eigenvectors 
\jm), where m = — j, — j + 1, . . . , + j. As usual we have 
the eigenvector relations j 2 \jm) = j(j + l)\jm) and 
jz\jm) = m\jm). Note that j can take on integer or half- 
integer vàlues, j = |, 1, |, . . ., such that N = 2, 3, 4, . . . 

The group of proper rotations in three dimension is 
denoted by SO(3) . This is the group of orthogonal 3x3 
matrices with determinant 1. An irreducible representa- 
tion of this group on the state space Tí is obtained in the 
Standard way: Given a rotation R e 5*0(3) the corre- 
sponding transformation of state vectors is provided by 
the unitary matrix 

D(R) = D(ni,n 2 , n 3 ) = exp (-in ■ jj . (2.1) 

The rotation R is characterized here by the vector n = 
(711,712,^3), i- e., R is the rotation about the axis given 
by n by the angle \n\ (in a right-handed sense). It should 
be mentioned that Eq. Ij2.ll) generally yields a two-valued 
representation of the rotation group: For half- integer j 
one obtains two unitary matrices which represent a given 
rotation R and which differ in sign. 



B. Rotational invariance of bipartite systems 

The representation Ij2.ll) leads to a representation of 
the rotation group on the tensor product space Ti (g> H of 
the bipartite system. If p is an operator acting on the ten- 
sor product space, a rotation R carried out on both parts 
of the composite system leads to the transformed oper- 
ator p' = [D(R) ® D(R)]p[D(R) (g> D(R)]ï, An operator 
p is said to be rotationally invariant or SO (3)-invariant 



3 



if it is invariant under all such transformation, that is, if 
the relation 

[D(R) ® D(R)]p[D(R) <g> D(R)]i = p (2.2) 

holds for allüe SO{3). 

A state of the bipartite system is given by a density 
matrix p satisfying p > and trp = 1. The set of all 
states p which fulfill the invariance requirement 112. 2|l will 
be denoted by S. It is clear that S is a convex subset of 
the set of all states of the bipartite system. 

The angular momentum operator of the composite sys- 
tem is given by J = j ® I + I ® j, where I denotes the 
unit matrix. The components of J are the generators 
of the product representation and the requirement of ro- 
tational invariance is equivalent to the statement that p 
commutes with all components of J . 

The product representation D(R) ® D(R) is obvi- 
ously reducible. Its decomposition into a sum of irre- 
ducible representations is a Standard subject of quan- 
tum mechanics. One introduces an orthonormal basis 
in Tt ® Tí which consists of the common eigenvectors 
\JM) of J 2 and J3 corresponding to the eigenvalues 
J(J + 1) and M, respectively, where J = 0,1,..., 2j 
andM = -J,-J+l,...,+J. The (2 J + l)-dimensional 
space which is spanned by the basis vectors \JM) with 
a fixed J is an invariant and irreducible subspace of the 
tensor product representation. 

The set S of rotationally invariant operators can now 
easily be characterized. To this end, we introduce pro- 
jection operators 



onto the angular momentum singlet state 



+J 



Pj= \JM){JM\ 



(2.3) 



M=-J 



which project onto the subspaces belonging to a fixed J. 
From the irreducibility of the representation within these 
subspaces one concludes with the help of Schur's lemma 
that any rotationally invariant operator p can be written 
as a linear combination of the projections: 
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(2.4) 



Here, the aj are c-numbers and we have introduced nor- 
malization factors (Ny/2J + It will be seen in 

Sec. IIII Dl that this choice of normalization factors leads 
to highly symmetric transformation properties of the pa- 
rameter space. For p to be Hermitian the aj must of 
course be real. Equation 12.411 then corresponds to the 
spectral decomposition of p. If p is a density matrix the 
aj are real and positive, aj > 0. On using trPj = 2J+1, 
the normalization condition takes the form 



2j 



j=o 



V2J+ 1 
N 



OíJ = 1. 



(2.5) 



For example, setting a — N and aj = for J = 
1,2, . .., 2j we get p = Pq = 1 00) (00 1 , i. e., the projection 



100) = -±= y (-1 

y/N ^ 



ú-m. 



li, m) ® |i, -m) 



(2.6) 



This state is the only pure state in S and it is maxi- 
mally entangled (the quantity ^ is the singlet fraction). 
Using the completeness of the projections Pj one con- 
cludes that the state corresponding to aj — \J2J + 1/N, 
J = 0,1,..., 2i, is the separable state p = j^I ® / of 
maximal entropy. 

It follows from the irreducibility of the representation 
D(R) that for any 5'0(3)-invariant state p the reduced 
density matrices p^ = tr 2 p and p^ = trxp of the sub- 
systems, given by the partial traces tr2 and tri, are pro- 
portional to the identity /. The reduced density matrices 
obtained from a rotationally invariant state thus describe 
states of maximal disorder. 

Summarizing, by means of Eq. 112.411 any rotationally 
invariant Hermitian operator is uniquely characterized by 
N real parameters aj. We can therefore identify the set 
of all such operators with the set of points 



a = 



(2.7) 



V Oí2j / 



in an 7V-dimensional parameter space R . The set of 
points a in this space satisfying aj > and the nor- 
malization condition H2.5H then describes the set S of 
rotationally invariant density matrices. In geometrical 
terms S represents an {N — l)-dimensional simplex. For 
instance, S is a line for N = 2, a triangle for N = 3, 
and a tetrahedron for N = 4. These examples will be 
discussed in Secs . IIVBl and WT% 



III. POSITIVE MAPS AND ROTATIONAL 
INVARIANCE 

A. Partial transposition 

Given an operator B on H the transposed operator 
TB = B T is defined in terms of the local basis states 
|Ím) by means of (jm\B T \jm') = (jm'\B\jm). Corre- 
spondingly, the partial transposition T2 = I ® T on the 
tensor product space is defined through 



T 2 (A (g> B) = A(g>TB = A® B 1 



(3.1) 



The operation of taking the partial transpose plays an 
important role in entanglement and quantum informa- 
tion theory. One reason for this fact is that T is a dis- 
tinguished example of a map which is positive but not 
completely positive E E El O El El • This means 
that T takes positive operators on H to positive opera- 
tors on Tí, while T 2 p need not be positive for a positive 
operator p on the tensor product space 7í®TL. 



4 



Important information on the entanglement structure 
of states is obtained by considering the action of positive 
but not completely positive maps. An example is given 
by the Peres PPT criterion according to which positivity 
under the partial tr ansp osition T 2 is a necessary condi- 
tion for separability An important general charac- 
terization has been developed by the Horodecki's [ÜJ]: A 
necessary and sufficient condition for a state p to be sepa- 
rable is that the operator {!®&)p is positive for any pos- 
itive map This condition does however not lead to a 
simple operational criterion for separability since we have 
no general structural characterization of positive maps, 
as it exists for completely positive maps in the form of 
the Kraus-Stinespring representation (lfil \WL lig] . 

B. $2-transformation 

If p is a rotationally invariant operator the partially 
transposed operator Tip is generally not invariant under 
rotations. It can be shown that, instead, T 2 p is invariant 
under transformations of the form D(R) (g) D(R)*, where 
D{R)* is the matrix obtained from D(R) by complex 
conjugation of its elements, that is D(R)* = D{R)^ T . 
Throughout this paper T denotes the transposition, f 
the adjoint and * the element-wise complex conjugation 
of a matrix. 

In the present investigation we shall utilize a map 
which is unitarily equivalent to the partial transposition, 
but which does map rotationally invariant operators to 
rotationally invariant operators. This map will be de- 
noted by i?2- By analogy to Eq. ll.'i.lll . $2 is taken to be 
of the form 

$ 2 {A®B) = A®ÜB = A®VB T V^ (3.2) 

with some fixed unitary matrix V. Hence, #2 = I ® $ is 
the partial transposition T 2 followed by a local unitary 
transformation acting on the second part of the bipartite 
system, that is, we have d 2 p = {I ®V)T 2 p{I ®V)^ . Since 
the maps $2 and T2 are unitarily equivalent a state p is 
obviously positive under d 2 if and only if it is positive 
under T 2 . 

The unitary matrix V will be determined from the 
condition that $2 preserves the rotational invariance of 
operators, i. e., if p is any invariant operator satisfying 
Eq. ll2.2ll we demand that the transformed operator ■& 2 p 
is again invariant: 

[D(R) <g> D{R)]d 2 p[D(R) <g> D(R)]Ï = ú 2 p. (3.3) 

This requirement is obviously satisfied if the map ú com- 
mutes with all rotations, that is, if the relation 

ti[D(R)BD(R)i] = D(R)(tfB)D(R)1 (3.4) 

holds true for all operators B on H and all R € 50(3). 
By use of the definition of $ given by Eq. Il.-S.2ll one can 
write Eq. <3.4ll as 

VD(R)*B T [VD(R)*]ï = D(R)VB T [D(R)V} 1 . (3.5) 



This equation is fulfilled if VD(R)* = D(R)V. Thus we 
see that the rotational invariance of -& 2 p follows from the 
rotational invariance of p provided we can find a fixed 
unitary matrix V such that 

VD(R)*V f = D(R) (3.6) 

for all R G 50(3). 

To obtain a unitary matrix V satisfying Eq. H3.6I) we 
employ specific properties of the representations of the 
rotation group. As in Eq. (12. H . let D(R) be the repre- 
sentation of the rotation R about an axis n = (n\, n 2 , ris) 
by an angle \n\. The complex conjugation of the elements 
of D(R) then yields the matrix 

D(R)* = exp(+ín-j T ^ 

= exp (—i[-niji + n 2 j 2 - mj 3 ]j 

= exp {-in' .j). (3.7) 

Here we use the fact that in the local basis \jm) the 
transposed components of the angular momentum op- 
erator are given by jf = j x , jj = -j 2 and = j 3 . 
Thus, D(R)* represents the rotation about the axis n' — 
(—ni, n 2 , —n 3 ) which is obtained from n through a rota- 
tion by 7r about the £2 -axis. To transform from D(R)* to 
D(R) we therefore define V to be the matrix represent- 
ing a 7r-rotation about the X2-axis. Using the notation 
introduced in Eq. <2.H we write 

y = D(0,7r,0). (3.8) 

Explicitly the matrix elements of V are given by 

(jm'\V\jm) = (-iy- m 5 m ,,- m . (3.9) 

Hence, V is real and we have V T = = V^ 1 . 
Equation l|3.8|l yields 

v(n'-j)V* = n.j, (3.10) 

which, by use of Eq. H3.7I) . immediately leads to the de- 
sired relation H3.6JI • We conclude that the map d 2 defined 
by Eqs. 13.211 and 13.8|l preserves the rotational invariance 
of operators. The advantage of this formulation is that 
i?2 , by contrast to T 2 , maps the set of rotationally invari- 
ant Hermitian operators onto itself and can be expressed 
as a simple linear transformation of the parameters aj. 
This transformation will be determined in Sec. IIII El 



C. Time reversal symmetry 

The transposition T is closely connected to the oper- 
ation of reversing the direction of motion, i. e., to the 
symmetry transformation of time reversal 0, O, l2ÏI |. 
We demonstrate that, in fact, it is the map d introduced 
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in Eq. )■'•>. 2t which describes the time reversal of partides 
with spin j. 

We have seen in the preceding subsection that T 
changes the sign of j% and leaves j\ and j'3 unchanged, 
while the unitary operator V (representing a 7r-rotation 
about the X2-axis) changes the signs of ji and j'3 and 
leaves j 2 unchanged. Hence, we have dj = Vj T V^ = —j. 
This shows that the map 1? describes the behaviour of the 
angular momentum operator under time reversal. 

It is known from Wigner's representation theorem 0] 
that the time reversal symmetry must be represented in 
terms of an anti-unitary operator. Indeed, we can express 
the action of 1? by means of an anti-unitary operator r 
through 



■dB = t-BV -1 . 



(3.11) 



The operator r = Vtq is composed of the unitary trans- 
formation V introduced above and of the anti-unitary 
transformation tq which is given by the complex conju- 
gation of the amplitudes in the basis \jm): 

\tp) = ^2c m \jm) 1 ^ r \ip) = ^2c* m \jm). (3.12) 

m m 

Thus, by virtue of Eq. A3.9J1 we have 

\<p) =Y,Cm\jm) » r\y) =Y,o* m {-l) j - m \j-m). 

rn m 

(3.13) 

This transformation expresses the well-known behaviour 
of spin-j partides under time reversal. For example, in 
the case N = 2 (j = ~ and m = ±|) Eq. IJ3.9I) leads to 
V = —i(T2, where a 2 is a Pauli matrix. The transforma- 
tion r thus consist of the complex conjugation and of the 
unitary transformation given by the matrix — Í02, which 
precisely corresponds to the time reversal transformation 
of a spin-i particle. 

In view of these results the map $2 = J ® 1? may be 
interpreted as a partial time reversal of the composite 
system. The fact that is not completely positive means 
that the operation of time reversal, when carried out only 
on a subsystem, does in general not lead to a physically 
legitimate state (22l |. 

D. Properties of the map $2 

The properties of the transformations $ and -d 2 are of 
course very similar to those of the transposition T and of 
the partial transposition T 2 , respectively. In particular, 
■d is a positive (but not completely positive) map, i. e., 
B > implies that i)B > 0. Moreover, -d preserves the 
trace, ti{i)B} = trB, and the unit matrix, dl = I. 

It follows from Eq. that 



(3.14) 



This equation illustrates the two-valuedness of the rep- 
resentation: V 2 = D(0, it, 0) 2 represents a rotation by 



2ir (i. e., the identity in SO(3)) and is equal to — / for 
half-integer j. Equation Ij3.14|l leads to the conclusion 
that, like T 2 , the map #2 is an involution which means 
that i?| — I ® ïï 2 ls equal to the identity map. In fact, 
employing Eq. H3.14JI we obtain for any operator B on tí: 

= V{VB T V J{ ) T V^ = VVBVÏVÏ = B. (3.15) 

Another property which will be important below is that 
i?2 is selfadjoint with respect to the Hilbert-Schmidt inner 
product, i. e., we have 



tï{x\d 2 Y)} =tr{(i? 2 ^)V} 



(3.16) 



for all operators X and Y on the tensor product space. 
This property derives from the corresponding property 
of the map 1?. Namely, for any two operators A and B 
on H we have according to Eq. H3.2J1 : 

tr{A\0B)} = tr{A^VB T V f } 

= tr{(V T A T V^ T ) J< B} 

= (-l) Aj tv{(VA T V^B} 

= tr{{iïA)ÏB}. (3.17) 

Note that we have used here that (— l) 4j = 1 for integer 
and half-integer j, and that V T = V^ 1 = (-l) 2j V, which 
follows from Eq. (I3.14jl . 

Since 1? is not completely positive the operator p' = í? 2 p 
need not be positive for a positive p. It is, however, 
invariant under rotations and can be represented in the 
form (12 . 4|) . To determine the action of í?2 on the aj- 
parameters we therefore write: 



1 23 

-Y 



a 



K 



1 23 



--Pi 



K 



a K 



d 2 P K , (3.18) 



where the parameters aK correspond to p and the a' K 
correspond to p' . We multiply this equation by Pj and 
take the trace using PjPk — SjkPk- This yields a linear 
transformation from the parameters aK to the parame- 
ters a' K . Using matrix notation we find 



6a, 



(3.19) 



where we have introduced a matrix with elements 



6 



1 



JK 



y/(2J+l)(2K+ï) 



tv{Pjiï 2 P K }. (3.20) 



The map z? 2 thus induces a linear transformation of the 
parameter space which is given by the N x N matrix 0. 

The matrix 6 is real symmetric, 6 T = <d, and orthog- 
onal, T 8 = /. The symmetry follows immediately from 
definition (13 . 20fi and the property ffi.lfijl . Since d 2 is an 
involution the matrix 9 must also be an involution, that 
is & 2 = I . Together with the symmetry of we therefore 
have T 8 = G 2 = /, which proves that is orthogonal. 
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E. Relation to Wigner's 6-j symbols linked to Wigner's 6-j symbols. To this end, we use 

Eq. I|3.20l) as well as the definition l|2.3f) of the projec- 
We derive a general expression for the elements of the tions P J in terms of the eigenbasis | JM), which gives 
matrix 0. It will be shown that these elements are closely 



e 



i 



+j 



-K 



JK 



= J2 E {JM\d 2 {\KQ)(KQ\)\JM). 

J- i n t t r~\ v 



(3.21) 



V(2J+l)(2íf+l) m —_jq—_jç 
To evaluate the $2-transformation in this expression we insert complete sets of product basis states \jmjm!) to get 

Qj K = —f= j== = E E (JM\iï 2 (\7n 1 m 2 )(m 1 m 2 \KQ)(KQ\m 4 m 5 )(m 4 m 5 \)\JM). 



y/(2J+l)(2K- 



M.Q m 1 ,m 2 m 4 ,m 5 



Here and in the following we shall frequently abbreviate \jmijm 2 ) by \mim 2 ), etc. According to the definition of $2 
[Eqs. and and to Eq. (EU we have 



l d 2 (\mim 2 ) (m 4 m 5 |) = \m\){m4,\ (g) V (\m 2 )(m 5 \) T = |mi)(m 4 | (gi V\m 5 )(m 2 \V' í 

= |mi)(m 4 | <g> (-l) 2 ^ m ^| - m 5 )(-m 2 | 
= (-i) 2 ^-"" 5 | mij _ m5 )( m4 , _ m2 | ; 



(3.22) 



which leads to 
Qjk 



vW+ïpTIj 



ÏÏ E E £(-u 



Af.Q mi,m 2 m 4 ,m 5 



2 ^ m2 "" l5 (JAf|m 1 ,-m 5 )(TO 1 m 2 |ií'Q)( J ft:g|m4TO 5 )(TO4,-m 2 |JM). 

(3.23) 

I 



The matrix elements in Eq. Q3.23I1 are vector-coupling 
(Clebsch-Gordan) coefEcients. Throughout the paper we 
adopt the usual phase conventions for these quantities, 
as they are given, e. g., in Ref. (Ï3 |. 

To evaluate further Eq. H3.23J1 it is convenient to em- 
ploy the 3-j symbols 



Ji J2 J3 
m\ m 2 7713 



(3.24) 



introduced by Wigner. These quantities are known from 
the theory of angular momentum coupling and are closely 
related to the vector-coupling coefEcients. Here, we have 



(mim 2 \JM) = (-1) M V2J+Ï 



3 3 
m\ m 2 



J 

-M 



(3.25) 

The 3-j symbols have many symmetry properties. The 
symmetry to be used here is given by 



3 3 J 
wi\ m 2 ?Ti3 



3 3 J 
-m\ —m 2 — 7713 

(3.26) 

We further need the selection rules for the 3-j symbols, 
namely that H3.24J1 is equal to zero for m 1 +m 2 + m 3 7^ 0. 

We introduce the relations (|3.25j) into Eq. (I3.23j) which 
yields a sum over products of four 3-j symbols. In the 



resulting expression we carry out the following manipula- 
tions: (i) we interchange the summation indices m 2 and 
TO5, (ii) we replace the summation index m\ by —mi, 
(iii) we introduce the new notation M = 1713, Q = mç, 
and (iv) we employ the symmetry relation l|3.26|) in the 
first and the third 3-j symbol. These manipulations lead 
to 

Q JK = V(2J+1)(2K+Í) J2 *( m *) ( 3 - 27 ) 

mi ,...,7716 

3 3 K 



x 



3 3 J 
m,\ m 2 m-i 

3 3 K 
— TO4 — m 2 m,Q 



m\ 777,5 —m§ 
3 3 J 

777-4 — 777-5 — 777-3 



where all sign factors have been collected in the quantity 



x(toí) = (-1) 



2 j — 771 2 — rn 5 + J + K 



(3.28) 



Finally, we use the selection rules for the first and the 
third 3-j symbol in Eq. Ij3.27l) which leads to m 1 +m 2 + 
n%3 = and —7714 — m 2 + m§ = 0. With the help of these 
relations it is easy to show that the phase factor x( m i) 
may be written as 

r-^y+mi ^_2) J+m2 (—]} J+ms 
x ^_ 1 y+m^_ 1 y+ m ^_ 1 ^ ) K+m e _ (3 39) 



x(m) 
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On using Eq. I}3.29H we see that the sum of the right- 
hand side of Eq. 13.27}! is exactly equal to a certain 6-j 
symbol of Wigner. The 6-j symbols are scalar quanti- 
ties which arise in the construction of invariants from 
the vector-coupling coefficients involving six angular mo- 
menta 0]. A general 6-j symbols is written as 



3l 32 33 
Í4 35 36 



(3.30) 



For the sum of Eq. 13.27}! we have j\ = j% = ji = j§ = j, 
ja = J and je — K . Hence, we finally obtain: 



Qjk = y/(2J+l)(2K+l) 



3 í J 
3 3 K 



(3.31) 



This equation represents a central result of this paper. It 
yields a general expression for the i^-transformation in 
terms of Wigner's 6-j symbols on which our investigation 
of the structure of rotationally invariant states is based. 

The 6-j symbols l}3.30ll are known to be invariant under 
any permutation of their columns and under the inter- 
change of the upper and lower entries in any two columns. 
It follows that the expression on the right-hand side of 
!}3.31}l is symmetric with respect to the interchange of 
J and K. It is also known from the theory of angular 
momentum that the expression on the right-hand side 
of 13.31}! represents an orthogonal matrix, in accordance 
with our previous considerations. 

The properties of the 6-j symbols have been studied 
in great detail and many explicit expressions and closed 
formulae are known. Computational methods and recur- 
sion relations for the 6-j symbols may be found in [ï^j . 
Equation l}3.31ü enables one to employ these results in 
the determination of the matrix O. For example, the 
first two rows and columns of 9 are given by 



N 



©JO 

Qji = Qu = y/3(2J + 1) 



(3.32) 

(JV-1)(JV + 1)-2J(J + 1) 



N(N- 1)(JV + 1) 
x (_l)2j+i+J. (3.33) 

Being real symmetric and orthogonal, the matrix O 
can of course be diagonalized and has eigenvalues ±1. 
The eigenvectors of may be found from the sum rules 
for the 6-j symbols given in ^4[. H we write the sum 
rule involving products of two 6-j symbols in terms of 
the matrix elements Qjk we get 



K 



ej K (-i) K e KL = (-i) L (-iy J e JL . 



(3.34) 



We infer from this equation that the vector oS L ^ with 
components a^' = (— 1) j 0jl is an eigenvector of with 
eigenvalue (— 1) L . Once we have determined the matrix 
we can therefore immediately write its eigenvectors: 
One multiplies for all J the J-th row of by (— 1) J ; 



the columns of the resulting matrix then represent the 
eigenvectors of 0. 

It follows from the orthogonality of the matrix that 
the vectors oS L \ L = 0,1,..., 2j, form an orthonor- 
mal basis of the parameter space. After a transfor- 
mation to principal axes therefore takes the form 
diag(+ 1 , — 1 , + 1 , . . . , ( — 1 ) 2j ; ) , which describes a refiection 
of the principal axes belonging to the eigenvalue —1. The 
trace of is obviously equal to zero for N even (half- 
integer j), and equal to 1 for N odd (integer j). 

According to Eq. !}3.32}l the components of the first 
eigenvector et' ' are given by af^ = (— 1) j 0jo = 
(-1) 2 · J V2J + Ï/N. This vector is proportional to the 
vector which represents the state of maximal entropy. 
The first eigenvector equation 9a' ' = a' ' thus ex- 
presses the invariance of the state of maximal entropy 
under It may be written as 



2j 



K=0 



JK- 



VïïkTÏ V2J + i 



N 



N 



(3.35) 



This equation can also be used to check that preserves 
the normalization J2.ófl . 



IV. SO(3)-INVARIANT PPT STATES 

A. Geomètric representation 

We define S p to be the set of 5'0(3)-invariant PPT 
states, i. e., the set of rotationally invariant states which 
are positive under i?2 (or, equivalently, under T^), The 
properties of $2 imply that S p is the set of density ma- 
trices p for which p' = -&2P is again a density matrix and 
that S p is the intersection of S with its image under #2 : 



s p = snti 2 s. 



(4.1) 



Since S is an (N — l)-simplex and $2 is a non-singular 
transformation the set 1D2S is again an (N — l)-simplex. 
Being the intersection of two convex sets, S p is also a 
convex set. 

With the help of the properties of the matrix derived 
in Sec. IIIIDl it is easy to give the geomètric construction 
of S p employing the space of the aj-parameters: One 
takes the (N — l)-simplex describing S and determines 
the intersection with its image under the linear map given 
by the matrix 0. Since S is convex it suffices to deter- 
mine the images of the extreme points of S in order to 
construct ^S. 

To facilitate the geomètric visualization we shall use 
in the following an (N — l)-dimensional parameter space: 
A Hermitian and rotationally invariant operator of trace 
1 is characterized uniquely by (N — 1) real parameters 
(ao,a 2 , ■ ■ ■ ,a 2 j-i). This means that we eliminate the 
parameter a 2 j by means of Eq. 12.5}! which expresses the 
condition of unit trace. The state space S can then be 
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identified with an (N— l)-simplex in 1 which is given 
by the conditions: 



^2 n aj ~ otQ,...,(X2j-i > 0. (4.2) 

j=o 

B. Examples 

We illustrate the geomètric construction of the set S p 
of PPT states for N = 2, 3 and 4. It will be seen that 
S p is isomorphic to an (N — l)-dimensional cube. The 
matrix elements Qjk can by determined with the help of 
Eqs. H3.32I) and 113. 33|) and by use of the general properties 
of 9 described in Sec. lïIÏDl 
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2 <g> 2 systems 



In the simplest case N = 2 the total system consists 
of two partides with spin i = \ (two qubits) . The to- 
tal angular momentum thus takes the vàlues J = 0, 1 
such that we can use a single parameter ao to describe 
a rotationally invariant Hermitian operator of unit trace. 
The inequalities (|4.2I) yield < ao < 2. The space of 
rotationally invariant density matrices is therefore given 
by the interval (1-simplex) S = [0,2]. The matrix G is 
found to be 



9 = 



1 



-1 V3 
%/3 1 



(4.3) 



which is obviously symmetric, orthogonal and of trace 
zero. The condition 112. 5J1 gives a\ — -^(2 — ao) which is 
used to eliminate a\ from the transformation a! = Qa. 
One finds that $2 maps the point ao = to a' — 1 and 
the point ao = 2 to a = — 1. This yields &2S = [—1, +1], 
and, hence, we get the set of PPT states: 



s p = snd 2 s = [0,1]. 



(4.4) 



We note that for the present case of two dimensions 
the rotational invariance is equivalent to the invariance 
under all product unitàries U®U. The states constructed 
above are therefore identical to the Werner states of 2 ® 2 
systems. 



2. 3 ® 3 systems 

For N = 3 (two qutrits) we have j = 1 and J = 0, 1, 2. 
We can therefore use two parameters (ao,ai) to char- 
acterize a Hermitian and rotationally invariant operator 
with trace 1. Equation 114.211 now yields that the set S of 
invariant states is given by the inequalities: 



1 



1 



7^0 



V3 



«i < 1, 



o;q, ai > 0. 



(4.5) 



Figure 1: SO (3) -invariant Hermitian operators of trace 1 for 
3®3 systems. Triangle ABC: The set S of rotationally invari- 
ant density matrices. Triangle A'B'C': The transform &2S. 
Rectangle ADA'E: The set S p of PPT stat es, wh ich is equal 
to the set S s of separable states (see Sec. IV C~2)l . The line 
DE represents the fixed points of $2, the point F the state of 
maximal entropy and B the singlet state |00). 



Hence, S is a triangle (2-simplex) with vertices A 
(0, 0), B = (3, 0) and C = (0, V3). 
The matrix 9 now becomes: 



9 = 



1 




(4.6) 



One easily verifies that this is a symmetric and orthogo- 
nal matrix of trace 1. On eliminating the parameter ct2 
we find that §2 acts as follows on the vertices of S : 



A = (0,0) 
B = (3,0) 
C = (0,V3) 




(4.7) 
(4.8) 
(4.9) 



Thus, ^25* is the triangle with vertices A' , B' and C , 

The sets S and #2 5 are depicted in Fig-Q] The figure 
also shows the line of the fixed points of #2 with end- 

points D = (1,0) and E = \ 0,^f)- This line is easily 
determined from the matrix 8 and its eigenvectors. Be- 
ing invariant under $2, the point F, which describes the 
state of maximal entropy, lies of course on this line. 

The rectangle with vertices A, D, A! and E represents 
the intersection S p — Sni^S 1 of the PPT states. It should 
be noted that the rotational invariance in the present 
example is equivalent to the invariance under the product 
transformations 0®0, where O varies over the group of 
orthogonal 3x3 matrices 0. 
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i 4 systems 



The case N — 4 corresponds to a system composed 
of two partides with spin j = |. The total angular 
momentum assumes the vàlues J — 0,1,2,3. Thus we 
get a three-dimensional parameter space with parameters 
(ao, a-y, «2). By virtue of Eq. I|4.2I) the set of rotationally 
invariant states is determined by the inequalities: 



1 V3 V5 

-a + —ai + —«2 < 1, 



ao, 0-1,0.2 > 0. (4.10) 



This shows that S is a tetrahedron (3-simplex) with ver 
tices A = (0,0,0), B = (4,0,0), C 



0,^5,0) and 



vertices of the tetrahedron 5 under $2: 

A= (0,0,0) " - 3 ^ 1 
B = (4,0,0) 



C 



D = 0,0, 




V5 



The points A' , B', C and D 1 are the vertices of the 
transformed tetrahedron &2S, as shown in Fig. EI 

We see from Fig.[21that the intersection S p = Síl^S is 
isomorphic to a 3-dimensional cube. An enlarged picture 
of this cube is shown in Fig. El The vertices of S p are 
given by the points A, A' and 



CM 




I^O.OJ. E' = 



G = 0,0, 



VE, 



G' = 



'2 2V3 



.0 



(4.16) 



(4.17) 



These points may be obtained as follows (see Fig.EJ. One 
takes the three edges emerging from the vèrtex A' of the 
tetrahedron Ú2S and determines their intersection with 
the faces of the tetrahedron S. This yields the points E' , 
F' and G' . The points E, F and G are then given by the 
images of E' , F' and G' under §2- 



Figure 2: S l O(3)-invariant Hermitian operators of trace 1 for 
4 ® 4 systems. The tetrahedron ABCD (continuous lines) 
represents the set S of invariant states, and the tetrahedron 
A'B'C'D' (broken lines) its transform -&2S. The intersection 
(bold lines) is the set S p of the PPT states. 



The matrix 6 is given by: 



= 



/ -1 

-VI 



Vï -VI V7 \ 



n 

5 

3 

5 

3V2T 



3V2Ï 
5 

7 

o 



(4.11) 



One checks that this matrix is symmetric, orthogonal and 
of trace zero. It leads to the following mapping of the 




1 



Figure 3: Enlarged picture of the cube S p of PPT states (see 
Fig-EJ. The plane AA'FF' subdivides S p into two prisms. 
The prism AA ' FF'G G' represents the set S s of separable 
states (see Sec. IV C~3ll . 
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V. SEPARABLE STATES 

A. Construction of S f O(3)-invariant separable states 

The set of separable states is defined to be the set of 
states p which can be written as a convex sum of product 
states: 



P : 



5> 



,(2) 



>0, = (5.1) 



where the p+ and pf^ are normalized local states Q]. It 
follows from this definition and from the positivity of "d 
that the map $2 is positive on separable states. Thus, 
i?2 maps rotationally invariant and separable states to 
rotationally invariant and separable states. 

We denote the set of S'0(3)-invariant separable states 
by S s ■ This set is contained in the set of states which are 
positive under #2: 



s s c s v = sn& 2 s. 



(5.2) 



This equation expresses the Peres PPT criterion. It can 
easily be applied in the present formulation once the ma- 
teix 6 has been determined: Given a rotationally invari- 
ant state p in terms of its parameter vector a according 
to Eq. 112.411 . a necessary condition for this state to be 
separable is that all components of the transformed pa- 
rameter vector a! = Qa are positive. 

To fully characterize the set of separable states one 
introduces a projection super-operator II, also known as 
twirl operator. Given any state p of the bipartite system 
the operator 



J=0 



(5.3) 



is positive, of trace 1 and rotationally invariant. The 
map p 1— ► IIp defines a projection (i. e., II 2 = II) from 
the total state space onto the space S of rotationally 
invariant states. Moreover, if p is separable then Up is 
again separable. 

We see from Eq. Ij5.3l) that the aj-parameters cor- 
responding to the projection Up are given by aj = 
-^===tr{Pjp}. If we take a pure product state 



p = I^ ( V 2) )V ( V 2) I 



(5.4) 



involving normalized local states l^ 1 -*) and the 
aj-parameters of its projection are found to be 



y J = àj[ í pM,<pW] = 



V2J+T 



v (^VV^ ( V 2) >. 

(5.5) 

It is known that any separable state can be written 
as a convex sum of pure product states. We define W 
to be the range of the parameter vector a whose compo- 
nents aj are given by the above functionals àj[(p^ , f^], 



where \<p^) and |</?*· 2 - ) ) run independently over all nor- 
malized states in H. With this definition one has the 
following result The set S s of rotationally invariant 
and separable states is equal to the convex hull of the 
range W, i. e., to the smallest convex set containing W, 
Thus, we have: 



S, = huU(W) C S p 



(5.6) 



The determination of S's therefore amounts to the deter- 
mination of the convex hull of the range of the functionals 
ttj^' 1 ',!/?' 2 )] given by Eq. l|5.5Jl . This task can be sim- 
plified by the following observations. 

First, since S s is the convex hull of W which, in turn, is 
contained in S p , a good starting point is to consider the 
extreme points (vertices) of S p , If one finds, for example, 
that all extreme points of S p belong to W one concludes 
immediately that S s must be identical to S p . 

Second, it is clear by construction that the functionals 
àj are invariant under simultaneous rotations l^ 1 ' 2 - 1 ) 1— > 
£>(i?)|<^ 1,2 )) of the input arguments. Pairs of state vec- 
tors differing by such a transformation are thus projected 
to one and the same point of the parameter space and 
need not be considered separately. 

Third, the range W is invariant under the map f} 2 . 
This means that if the point a belongs to W, then also 
the transformed point Qa belongs to W. This statement 
can easily be proven by use of the results of Sec. IIII Cl 
In fact, we have 



àj[<pW,T<pW] 



N 



\/2J+ 1 

2j 



(^V 2) I^jMV 2) ) 

(5.7) 



K=0 

We see that the transformation a. 1— * Qa corresponds 
to the time reversal transformation r carried out on 
the second input argument of the functionals. Equation 
115. 7|l also demonstrates that if |<^ 2 )) is invariant under 
r the corresponding parameter vector represents, for any 
choice of l^ 1 )), a fixed point of G. 



B. Representation in terms of spherical tensors 

In addition to the projections Pj there exist further ro- 
tationally invariant operators which span the set S and 
which lead to a particularly useful representation of the 
set of separable states. To construct these operators 
we introduce the irreducible spherical tensor operators 
Tjm acting on Tl, where, as before, J = 0, 1, . . . , 2j and 
M = — J, — J + 1, . . . , + J. The matrix elements of these 
operators are defined by the 3-j symbols: 

3 í J 



(jm\T JM \jm') = (-iy- m V2J + ï 



m — m —M 



(5.8) 

According to the selection rules of the 3-j symbols the 
matrix element 115. 8JI is zero for Am = m — m' ^ M . 
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The tensor operators Tjm represent a complete sys- 
tem of operators on Tí which are orthonormal with re- 
spect to the Hilbert-Schmidt inner product, i. e., one has 

to{TjMTj'M'} = àj.J'ÒMM'- 

For a fixed J the (2J + 1) operators Tjm transform 
according to an irreducible representation of the rotation 
group corresponding to the angular momentum J. For 
example, the T\m transform as the spherical components 
of a vector, while the T 2 m behave as the components of 
a second-rank tensor under rotations. For N = 2 the 
tensor components T\m may be expressed in terms of the 
Pauli matrices as T 10 = and ïi f ±i = t|(<ti ±ia 2 ). 

The definition l|5.8l) leads to the relation T JM = Tj M — 
{—l) M Tj_ M . One concludes that the tensor operators 
are eigen-operators of the time reversal transformation: 
§T JM = {-1) J T JM . 

It follows from the transformation behaviour of the 
Tjm that the operators on the product space defined by 



Qj= t jm®t] 



JM 



(5.9) 



where — 1 < A < +1. These states are invariant under 
all product unitàries U ®U. Therefore, all states of the 
family are, in particular, invariant under rotations and 

are found to be 



belong to S. The parameters corresponding to pw 



aj 



N 



V2J+ 1 
V2J+ 1 

Jsr* 



tr{P JPw } (5.15) 
X+(-l) 2j+J (NX-l)] . 



To obtain this result one has to determine the expression 
tr{P,/F}. This may be done by noting that for J = 
Eq. {Olt yields F = N<d 2 Po- The expression tr{P/F} 
can therefore be written in terms of the matrix elements 
6jo which are given by Eq. I|3.32l) . The family of the 
isotropic states can be embedded in a similar way into 
S if one first performs the local unitary transformation 
I®V. 



C. Examples 



are invariant under rotations. The connection between 
the projections Pj and the operators Qj is provided by 
the relation 



$2Pj = QjF, 



(5.10) 



We construct the set S s of separable states for the ex- 
amples considered in Sec. lïVBl To this end, we make use 
of the functionals H5.12J1 which characterize S s and of the 
general properties of the range W described in Sec. IV Al 



where we have introduced the flip operator F which is 
defined by 



¥\jmijm 2 ) = \jm 2 jm 1 



(5.11) 



The proof of Eq. H5.10I) is given in AnnendixlAl 

Equation Ij5.1()|l leads to an alternative characteriza- 
tion of the set of separable states. Since \<fi^} and Iv^ 2 ') 
vary independently over all normalized states we may 
use the right-hand side of Eq. I|5.7|) instead of the origi- 
nal expression 15. ót for the functionals àj[tp^\ f^}. If 
we introduce l|5.1()|) into ljB.7ll we find that we can employ 
the functionals 



aj[^\^ 



N 



VïïT 



(2)\|2 



M=-J 



(5.12) 

in order to construct the range W and the set S s of sep- 
arable states. An advantage of this formulation is that 
it leads to a very simple expression for J = 0. Namely, 
since Toq = -L·l we have 

V JV 



(5.13) 



It might be interesting to note that Eq. H5,10J1 can be 
used to identify the one-parameter family of the Werner 
states given by 



1 



Pw 



N 3 — N 



[(N - X)I <8> I + (NX - 1)F] , (5.14) 



i 2 systems 



In the case of our first example discussed in Sec. lIVBTl 
we found that the parameter ao describes a PPT state 
if and only if ao € S p — [0, 1]. We immediately see from 
Eq. I|5.13|l that the set of separable states and the set of 
PPT states are identical, that is S s = S p . In fact, accord- 



ing to Eq. H5.13Ü the functional ào^ 1 ^ 



(2)1 



can take any 



value in the interval [0,1] because |</?' 1:2 ^) are arbitrary 
normalized states. This shows that in the present case 
positivity under $2 is a necessary and suficient condition 
for separability, which is a well-known fact [23| . 



2. 3 ® 3 systems 



Using the results of Sec. HVB"2l we show that also for 
3 ® 3 systems the set of PPT states and the set of sep- 
arable states coincide, i. e., S's = S p . Thus, positivity 
under d 2 is again a necessary and sufficient condition for 
separability in this case, as has been demonstrated by 
Vollbrecht and Werner To prové this we verify that 
the extreme points of S p , that is, the points A, A' , D and 
E belong to the range W (see Fig.nj. 

First, we choose \ = \j = l,m= +1) and \<p^) — 
\j = 1, m = — 1). These states are orthogonal and, hence, 
à = according to Eq. H5.13J1 . Using the selection rules 
for the matrix elements H5.8Ü of the tensor operators T ÍM 
one sees that also &\ = (the operators T\m cannot 
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connect states whose magnètic quantum numbers differ 
by 2). This shows that A — (0,0) belongs to the range 
W. It also follows that A' belongs to W, because A' is 
the image of A under "d 2 • 
Next, we consider the state 

b (2) > = ^=(|l,+l) + |l,-l». (5.16) 

This state is invariant under the time reversal transfor- 
mation ()3.13f . Thus, for any choice of the point 

(ào, ài) is a fixed point of d 2 and, hence, belongs to the 
line DE (see Fig. QJ. If If^) is any state orthogonal 
to l^*- 2 )) we have that, additionally, ào = and, hence, 

(à ,ài) = (o, = E £ W. One the other hand, 

if we take |</5^) = \ Í P^ 2 '), then à — 1 and, therefore, 
(ào, ài) = (1,0) = D £ W . This concludes the proof. 



3. 4 ® 4 systems 

For A® 4 systems it is again possible to give a complete 
geomètric construction of the set of separable states by 
use of the results of Sec. IIV B 31 In the case N = 4 we 
have to consider the following functionals: 



à [^V 2) ] = KW )!'. 



h[<P (1 



= 71 £ 



(5.17) 
^M 2 ,(5.18) 



M=-l 
+ 2 



à2 [^),^)] = ^=J2 l(^ (1) |r 2 M|^ 2) )| 2 .(5.19) 

^ 5 M=-2 

To construct S s we proceed in four steps, investigating 
the extreme points of S p given by Eqs. (14. 1 2t . l|4.1fill . 
(ETTTjl and (gjg) (see Fig.®. 

(1) We show that A,A' £ W. To proof this we take 
\íç( 1 )) = ||,+|) and |</? (2) ) = ||,-|). These sta tes ar e 
orthogonal and, therefore, ào = accordingto Eq. Q5.17J1 . 
Since the magnètic quantum numbers of the states differ 
by 3 the selection rules for the matrix elements 115. 811 yield 
that ((p {1) \TjM\tp {2) ) = for J = 1,2. Thus, Eqs. (IB~T3j) 
and l|5.19ll yield àx = ct 2 = 0. Hence, A — (0,0,0) and 
A' = tf 2 A belong to W. 

(2) We demonstrate that also G, G' £ W. To this end, 
we take \<pW) = ||,+|) and |^ 2 )) = | § , — i> . These 
states are again orthogonal and we get ào = 0. Since 
Am = 2 the matrix elements (y^ 1 ) \Tim\^ 2 ^) vanish and, 
therefore, ài = 0. The only matrix element of the T 2 m 
which is not equal to zero on account of the selection 
rules is given by 



22W 



(2)\ 



1 

V2- 



(5.20) 



Thus, with the help of Eq. H5.19J1 we obtain à 2 = Sg. 
One concludes that G = (0, 0, -7=) and, hence, also G' = 
■d 2 G belong to W. 



(3) We claim that F, F' £ W. To prové this we choose 
the states: 



,(2)\ 



V2 ( 

= 7l( 



(5.21) 
(5.22) 



These states are obviously orthogonal and we get again 
ào = 0. The selection rules now yield (ip^ \Tj^i\(p^) = 
for J = 1, 2 and M ^ 0, while 

1 /3 3 



(2)\ _ 



2 \2' 
1/3 
2 \2 : 



T 



JO 



Tj 



|,-Í).(iu.) 



We have the following general relation between the ma- 
trix elements of the tensor operators: 

(j,-m\Tj \j,-m) = (-l) J (j,+m\Tj \j,+m). (5.24) 

On using this we see that the expression Q5.231 vanishes 
for J = 2. It follows that à 2 = 0. On the other hand, for 
J = 1 we obtain 



(^\T W \^ 2) ) 



Tu 



2^5' 



(5.25) 



With the help of H5.18jl this leads to ài = ^p. In sum- 

mary, we see that F = (0, ^0-, 0) and F' = -d 2 F belong 
to the range W . 

(4) It is shown in Appendix El that the functionals 
ll5~T7jl and fTTflt fulfill the inequality: 



à 2 [íp 



(5.26) 



It follows that E and E' do not belong to the range W, 
Namely, for these points we must have a 2 — and ào = | 
[see Eq. H4.16J1 ] which contradicts the inequality l|5.26|l . 
This shows that in the present case S's is a true subset 
of S p , i. e., positivity under d 2 is a necessary but not 
sufücient condition for separability. 

According to Sec. IIV B 31 the set S's of separable states 
is contained in the cube S p of the PPT states (see Fig.EJ. 
By the above results the points A, A' , F, F', G, and G' 
are contained in the range W. Since S s is the convex hull 
of W we conclude that S s contains at least the polyhe- 
dron AA' FF'GG' . We observe that this polyhedron is 
isomorphic to a prism. 

The inequality H5.26H yields an additional condition for 
the separable states. It implies that all points of the 
range W must lie on or above the plane which is defined 
by a 2 = -^«0 and which is indicated as gray surface in 
Fig.0 We note that according to Eqs. I|4.12|) and 14,17|) 
the points A, A' , F and F' belong to this plane. It follows 
that the set S s of separable states is in fact identical to 
the prism AA' FF'GG' . 
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In summary, the convex structure of the set of SO{3)- 
invariant states of 4 ® 4 systems may be described by the 
following inclusions: 

(prism S s ) C (cube S p ) C (tetrahedron S). (5-27) 

The tetrahedron S, representing the set of all invariant 
states, decomposes into the cube S p of PPT states and 
the set S \ S p of entangled states whose partial trans- 
position has negative eigenvalues. The cube S p of PPT 
states, in turn, consists of the prism S s of separable states 
and of the set S P \S S of entangled PPT states. The plane 
Q'2 = 4g ao thus separates the entangled PPT states from 
the separable states. 

As can be seen from Fig.Elthe set S P \S S is isomorphic 
to a prism from which one face has been removed. All 
states belonging to this set are inseparable and have pos- 
itive partial transposition. This leads to the important 
conclusion that S p \ S s represents a three-dimensional 
manifold of bound entangled states, i. e., states which 
cannot be distilled by loc al q uantum operations and clas- 
sical communication [ÏTllÜ 25 1. 



VI. DISCUSSION AND CONCLUSIONS 

We have analyzed the structure of the state spaces of 
bipartite N®N systems which are invariant under prod- 
uct representations of the rotation group. The main tool 
of the analysis is the positive map d which is unitarily 
equivalent to the transposition T and describes the be- 
haviour of local states under time reversal. Employing 
the properties of $ one relates the partial time reversal 
i?2 = I <8> i? to a linear transformation of the parameter 
space = {ex} and expresses the corresponding matrix 
6 in terms of Wigner's 6-j symbols. This matrix has 
been used to obtain geometrical representations for the 
sets of the separable and of the PPT states in the cases 
N = 2, 3 and 4 

In Sec. IV C 31 the inequality H5.26J1 enabled the con- 
struction of the set of separable states. Taken together 
with the Peres PPT criterion this inequality yields a nec- 
essary and sufficient condition for the separability of ro- 
tationally invariant states of 4 <g> 4 systems. It is of great 
interest to examine the possibility of an extension of this 
picture to higher dimensions. In this context it is impor- 
tant to observe that the inequality ljFi.261) expresses the 
positivity of a certain map $ which is given by 



<S>B = 



+2 

£ 

M=-2 



T 2 mBT 2M — T oa BT OQ . 



(6.1) 



This map is non-decomposable and deteets all entangled 
PPT states. Hence, we need exactly two maps, namely 
and <£>, in order to identify uniquely all separable states. 
These maps yield complementary conditions for separa- 
bility in the sense that the two inequalities 



constitute a necessary and sufficient separability crite- 
rion. It should also be noted that the proof of Appendix 
[5]does not rely on any invariance requirement. We con- 
clude that positivity under the map $2 = 7 <8> $ is a 
necessary condition of separability for all (not necessar- 
ily rotationally invariant) states of 4 ® 4 systems. 

The positive map introduced in Eq. Hfi.ll) corresponds 
to an entanglement witness |2(1 12(1 l27l | which is given 
by the operator W = P2 — Pq. The plane «2 = "^f a o 
in parameter space may be viewed as an optimal hyper- 
plane defined by this witness W. This fact leads to the 
following interpretation of the inequality (|5.2fi|) : If a mea- 
surement of the total angular momentum is carried out 
on a separable state, the probability of finding the value 
J = 2 must be larger or equal to the probability of finding 
the value J — 0. 

The method developed here suggests many general- 
izations and applications. An obvious extension is to 
consider bipartite systems whose local state spaces are 
not isomorphic, involving two different angular momenta 
jW ^ j( 2 ). Further important tòpics are an extension 
of the analysis given in Sec. IV Cl to higher-dimensional 
systems, the treatment of other symmetry groups, and 
entanglement in multipartite systems. 

The matrix O contains the complete information on 
the behaviour of the spectrum of the invariant states un- 
der partial transposition. It can also be used to express 
various separability criteria and entanglement measures 
and to design positive maps and entanglement witnesses. 
Examples of applications are the determination of the 
relative entropy of entanglement with respect to the set 
of PPT states [2fl] an d the entanglement measure given 
by the negativity |2JJ, [3JJ|. The negativity, for instance, 
is determined by the trace norm of the partially time- 
reversed state which can be written as 



>p|ii = E 



V2J+1 

N 



K 



(6.3) 



{I®ú)p> 



and 



{I®$)p > 



(6.2) 



where \\A\U = tr\A\ denotes the trace norm of A. 

In Refs. y| Lü| a necessary separability criterion, the 
reduction criterion, has been introduced which is based 
on the positive map defined by AB = ItïB — B. This 
criterion is not stronger than the Peres criterion, but has 
the important benefit that any state violating it can be 
distilled. For SO (3)-invariant states the reduction crite- 
rion is equivalent to the inequality based on the quan- 
tum Renyi entropy Soa S, Üj] and to the disorder cri- 
terion and takes the form -^7 ® I — p > 0. In 
terms o f the pa rameters a,j this can be expressed through 
atj < V2 J + 1. We see explicitly from our examples that 
for rotationally invariant states the reduction criterion 
is in fact much weaker than the Peres criterion. For in- 
stance, in the case N = 4 we get from it the conditions 
«0 < 1 and ai < v3. The region defined by these in- 
equalities is much larger than S p and than the true set 
S s of separable states (see Fig.|5||. 

Recently, a necessary criterion for separability has been 
developed by Rudolph [33I l34| , which is known as cross 
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norm or realignment criterion [35] . This criterion is 
based on the cross norm of the states of the tensor prod- 
uct space flfl and provides strong conditions for separa- 
bility. It is generally neither weaker nor stronger than the 
PPT criterion. It can detect, however, bound entangle- 
ment. To formulate the cross norm criterion we associate 
with any density matrix p — J2i Q ® D i a map $ p by 
means of the formula 



$ j9 S = ^C i tr{(í?A)S}. 



(6.4) 



For a separable state p the corresponding map $p is a 
contraction with respect to the trace norm, i. e., we have 
H^plli _• 1) which immediately yields a necessary condi- 
tion for separability. 

The application of the cross norm criterion to rota- 
tionally invariant states leads to an inequality which can 
again be expressed entirely in terms of the elements of the 
matrix 8. If the state p is given by its parameters aj 
the trace norm of $ p can be written in a form analogous 
to Eq. 



1$, 



E 



V2J+ 1 
N 



A) K a K 



K 



< 1. (6.5) 



This is a general expression for the cross norm criterion of 
50(3)-invariant states in any dimension N. It allows an 
explicit determination of the regions in parameter space 
satisfying or violating the criterion. In particular, with 
the help of the above formula one immediately evaluates 
the trace norm | |$ p | |i for the families of the Werner states 
and of the isotropic states. 

We finally mention that the present results could also 
find a number of important applications in the theory of 
open systems j^. The close connection to open system 
is based on an isomorphism [38] between states p on the 
tensor product space TL®TL and completely positive maps 
$ of operators on 7i. We define this isomorphism by the 
relation p = (7®$)Po. Apart from a normalization factor 
this relation is equivalent to Eq. llíi.4ll . It yields a one- 
to-one correspondence between the rotationally invari- 
ant density matrices p and the completely positive maps 
$ which are trace-preserving and rotationally invariant. 
Such maps arise through the interaction of open systems 
with isotropic environments. The isomorphism thus al- 
lows one to use the structure of S in the construction of 
appropriate representations of one-parameter families of 
quantum dynamical maps and to derive the general form 
of isotropic non-Markovian quantum processes. 



Appendix A: PROOF OF RELATION lprïo|l 

In the basis of the product states \m1m2) = \jm1jm2) 
the matrix elements of the operator faPj are found to 
be 

(m 1 m 2 1 $2 Pj | rrí-y m' 2 ) 

= (-l) 2j -" l2 ""^(m 1 ,-m 2 |P J |m' 1 ,-m 2 ), (Al) 



(m,-m'\JM), (A2) 
[m',-m\JM). (A3) 



where we have used the definition 13.21 of the $2- 
transformation as well as the matrix elements 113. 9J) of 
the unitary matrix V introduced in Eq. 113. 8J1 . On the 
other hand, the definition l|5.8H of the tensor operators 
and Eq. Ij3.25j) lead to 

(m\T JM \m') = (-iy- m ' 
{m\T] M \m') = (-I)*'™ 

We recali that the matrix elements on the right-hand 
sides are vector-coupling coefficients which are taken to 
be real following the usual phase conventions. The defi- 
nitions H5.9J1 and H5.11H for the operators Qj and for the 
flip operator F yield: 

(mim 2 \Qj¥\m' 1 m' 2 ) 
+.J 

= E {-^f i ~ m2 ~ m ' 2 {rni,-m' 2 \JM){JM\m' 1 ,-m 2 ) 
M=-J 

= (-l) 2 3- m2 ~ m z(m 1 ,-m 2 \P J \m' 1 ,-m 2 ). (A4) 

Comparing this with l|Alj) we see that QjW — 1Í2P./, as 
claimed. 



Appendix B: PROOF OF INEQUALITY l|06)l 

We take any fixed normalized state \<p^) and decom- 
pose it with respect to the basis states |m) = \jm): 



k (2) )= Cl 

















+ C 2 




+ C 3 


-D 















(Bl) 



The normalization condition for the amplitudes Ci reads 

4 

En 2 = 1 - ( B2 ) 



Consider then the operator: 



+2 



A =^ E T2m\^ 2 W 2) \T, 



2M- 



(B3) 



M=-2 



This operator is obviously Hermitian and positive and 
we have à 2 [^ 1 ',^ 2 '] — (<p^ \A\tp^). It will be demon- 
strated below that l^*- 2 -*) is an eigenvector of A corre- 
sponding to the eigenvalue ^= : 



1 



This equation implies that A can be written as 
.4 



i+-L|^))(^)|, 



(B4) 



(B5) 



where A is again a positive operator. This leads to 

«2 



5 2 [^ (2) ] = (^|i|^) + -J=|(^|^)| 2 



> -l^V 2 )] 



(B6) 
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which proves the inequality 115, 26|) , tation of the operator A in the basis |m). With the help 

It remains to demonstrate the eigenvector relation of the matrix elements Ij5.8l) of the tensor operators Tjm 
l|B4|) . To this end, we determine the matrix represen- one flnds that A is represented by the matrix 



V5 



| Cl | 2 + 2|c 2 | 2 + 2|c 3 | 2 

-C\C2 + 2CgC4 

— c[cz — 2C2C4 



-C1C2 + 2c 3 cl 
| C2 | 2 + 2| Cl | 2 + 2|c 4 | 2 
C2C3 
— C2C4 — 2c|c3 



C3 



-C1C3 - 2C2C4 
C2C3 

+ 2|c 4 | 2 + 2| Cl | 2 
-C3C4 + 2c|c 2 



C4 



CiC 4 

-C2C4 — 2C1C3 

-C3C| + 2ciC2 

+ 2|c 3 | 2 + 2| C2 | 2 



(B7) 



It is now easy to verify by an explicit calculation that 
the vector c = [c\, c 2 , c 3 , C4) 71 , which represents the state 
|<^ 2 )) according to Eq. lIBlIl . is an eigenvector of this ma- 



trix corresponding to the eigenvalue -j=. This concludes 
the proof. 
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